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Abstract 

The saturated dynamics of a Single Pass Free Electron Laser is considered within 
a simplified mean field approach. A method is proposed to increase the size of the 
macro-particle, which is responsible for the oscillations of the intensity of the wave. 
This approach is based on the reconstruction of invariant tori of the dynamics of 
test particles. To this aim a dedicated control term is derived, the latter acting 
as a small apt perturbation of the system dynamics. Implications of these findings 
are discussed in relation to the optimization of the laser source. PACS : 41.60.Cr, 
47.10.Df, 05.45.Gg 
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Introduction 



A Free Electron Laser (FEL) generates a tunable, coherent, high power radi- 
ation. A FEL differs from conventional lasers in using a relativistic electron 
beam as its lasing medium. The physical mechanism responsible of the light 
emission and amplification is the interaction between a relativistic electron 
beam and a magnetostatic periodic field generated in the undulator. Due to 



1 Tel. : +33 4 91 26 95 47 ; Fax : +33 4 91 26 95 53 
Email : bachelard@cpt.univ-mrs.fr 

2 Unite Mixte de Recherche (UMR 6207) du CNRS, et des universites Aix-Marseille 
I, Aix-Marseille II et du Sud Toulon- Var. Laboratoire affilie a la FRUMAM (FR 
2291). 



Preprint submitted to Elsevier 



February 8, 2008 



the effect of the magnetic field, the electrons are forced to follow sinusoidal 
trajectories, thus emitting synchrotron radiation. This initial seed, termed 
spontaneous emission, is then amplified along the undulator until the laser 
effect is reached. 

Among different schemes, Single-Pass high-gain FELs are particularly attrac- 
tive since they hold the promise of resulting in flexible light sources of broad 
theoretical and applied interests. The coupled evolution of the radiation field 
and the N particle beam in a Single-Pass FEL can be successfully modeled 
within the framework of a simplified Hamiltonian picture [Tj . The N + 1 de- 
gree of freedom Hamiltonian displays a kinetic contribution, associated with 
the particles, and a potential term accounting for the self-consistent coupling 
between the particles and the field. Hence, direct inter-particles interactions 
are neglected, even though an effective coupling is indirectly provided because 
of the interaction with the wave. 

The theory predicts a linear exponential instability and a late oscillating sat- 
uration, for the amplitude of the radiation field. Inspection of the asymptotic 
phase-space suggests that a bunch of particles gets trapped in the resonance 
and forms a clump that evolves as a single macro-particle localized in space. 
The remaining particles are almost uniformly distributed between two oscil- 
lating boundaries, and populate the so called chaotic sea. 

This observation allowed to derive a simplified Hamiltonian model to charac- 
terize the saturated evolution of the laser. Such reduced formulation consists 
in only four degrees of freedom, namely the wave, the macro-particle and the 
two boundaries, delimiting the portion of space occupied by the the uniform 
halo surrounding the inner bulk [2][3] . 

Furthermore, the macro-particle rotates around a well defined fixed point and 
this microscopic dynamics is shown to be responsible for the macroscopic os- 
cillations observed at the intensity level. It can be therefore hypothesized that 
a significant reduction in the intensity fluctuations can be gained by imple- 
menting a dedicated control strategy, aimed at confining the macro-particle 
in space. In addition, the size of the macro-particle is directly related to the 
bunching parameter, a quantity of paramount importance in FEL context [3|. 

Smith et al [4J showed that a test-wave, being characterized by a frequency 
close to the one of the wave, can effectively destroy the macro-particle and 
activate a consequent detrapping process. Dimonte et al [5] implemented this 
approach on a Travelling Wave Tube and detected a significant reduction of the 
intensity oscillations, followed, however, by an undesired systematic collapse 
of its mean- value. 

In this paper, we focus on the macro-particle, spontaneously established as a 
result of the wave-particle interaction process in the saturated regime. In par- 
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Figure 1. Left : Normalized intensity of the FEL's radiation simulated from Hamil- 
tonian (HJ. Right : N = 10000 electrons in phase space, at t = 400, 401.25, 402.5 
and 403.75. 

ticular we develop a dedicated technique to influence and, possibly, control its 
evolution, thus opening up the perspective of defining innovative approaches 
aiming at stabilizing the laser signal. To this end, we consider a test-particle 
Hamiltonian in a mean-field approach and calculate a small but appropri- 
ate control term which acts as a perturbation. The latter is shown to induce 
an increase in size of the macro-particle. A regularization of the dynamics is 
also observed, as confirmed by the reconstruction of invariant tori around the 
massive core. 

The paper is organized as follows. In Sec. [fl we introduce the mean-field model 
and shortly outline its derivation from the original iV-body Hamiltonian. In 
Sec. El a test-particle is "controlled" , through the reconstruction of invariant 
tori of its dynamics. Finally, we draw our conclusion and discuss possible 
implication of the present analysis. 



1 The mean-field model 



As previously anticipated, the dynamics of a Single-Pass FEL is successfully 
captured by the following Hamiltonian [1] : 

^ = £[f -2^008(^ + 0)], (1) 

where (0i,Pi) are the position and momentum of the 2-th particle and (4>,I) 
stand respectively for the phase and intensity of the radiation. In the mean- 
field model, the conjugated variables (0, J) are replaced by two functions of 
time and I{t), obtained from direct simultations of the self-consistent 
dynamics. In other words, this amounts to formally neglecting the action of 
the electrons on the field. 
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The iV-body Hamiltonian ([T]) can therefore be mapped into 

N 

t=i 

where 

H lp (9,p, *) = y - 2 ^cos (6 + 0(t)). (2) 

Hence, the dynamics of a FEL can be adressed by studying the evolution of 
a test particle, obeying Hamiltonian (j2J). The functions I(t) and <p(t) act as 
external fields and are here imposed by assuming their simplified asymptotic 
behaviour as obtained by a frequency analysis [6] : 

k=l 

in the reference frame of the wave. 



2 Hamiltonian control of a test-particle 



The Hamiltonian control addresses systems which are close to integrable, i.e. 
whose Hamiltonian can be written as H = H + V, where H is integrable and 
V a perturbation of order e (compared to H ). The results we use here have 
been proven rigorously [7118] . In practice, it can be shown that a suitable control 
term f of order e 2 exists such that Ho + V + f has an invariant torus at a given 
frequency u . In our case, the perturbation corresponds to the oscillating part 
of the intensity. The interaction term of Hamiltonian (T5]) reads 



2 V C ° S ^ ~ ^ = FcosO — eRe(e i9 W(t)). 
Thus, our integrable Hamiltonian can be cast in the form 

2 

H = P --F cos9, (3) 

whereas the time-dependent perturbation V is identified as 

V{t,9) = eRe(e ie W(t)). (4) 

First, we express Hamiltonian ([3]) into action-angle variables (ip, J) [9]. Then, 
we expand Ho around J = Jo, which in turn identifies the region where the 
invariant torus is reconstructed in : 

H Q (J) = E + uj (J - Jo) + S(J - Jo) 2 + 0((J - Jo) 3 ). 
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Likewise, the ^-component of perturbation is expanded as 

M L 

e"=£ £ «m,n(J- J )"V^ + O((J- J ) M+1 ) 

m=0 n=—L 

The control term reads [8] 

/(</?, t) = V(J , t) - V{ Jo - Td v V(J , <p, t), <p, t), 
where T is a linear operator acting on an element of the Fourier basis as : 

p i(u)t+n<p) 



i{oj + nuo) 

At the second order in e (which is the main term), the control term reads 

f(cp, t) = e 2 w(v, t)Td v v(cp, t) - <?8{Td v v{ip, t)f, 
where v and w are the first terms in the expansion of V : 

V(J, <p, t) = v(<p, t) + (J- J )w(<p, t) + +0((J - Jo) 2 ), 

which gives 

k=0n=-L Uon + U k 

and ^ 

™fo*) = £ E 

fc=0n=-L 

where are the Fourier coefficients of : 

K 

W{t) = W k e iuJkt . 

Numerically, two Fourier modes are considered for W(t) (K = 2), and eleven 
for the ^-dependent part of the perturbation (J4j) (L = 5). The perturbation is 
of amplitude e = 1/5. The expansion is performed around Jo ~ 1-33 (which 
corresponds to the energy E — 0). 

In action-angle variables, the regularization of the dynamics (see Figj2]) is 
clearly seen. Rather than a single torus, the control term generates a contin- 
uous set of invariant tori, thus expanding the regular domain in phase space. 

Unfortunately, the exact change of variables from (if, J) to (9, p) [9] has a 
singularity at the pendulum separatrices. In order to implement our control 
on the whole space, we use a simplified, but regular, change of variables which 
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Figure 2. Poincare sections of a test-particle of Hamiltonian Hq(J) + V(J, tp, t) (left) 
and Hq(J) + V(J, cp, t) + f((p, t) (right), in the pendulum action-angle variables. 
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Figure 3. Left : control term /(</>, t), at t = 0. Right : comparison of the exact and 
simplified changes of variables : (p is plotted as a function of 9, for a given energy 
(and so for a given J). The plain line corresponds to the exact change, the dotted 
to the simplified one. The curves correspond, from bottom left to bottom right, to 
E = 0.5, and -0.5. 

mimics the exact one in the region of the invariant torus predicted by the 
control (the error in this region is less than 4%) (see Figj3]) : 

6 

tan ip = -. 
V 



In terms of the (9,p) variables, our control term now reads : 

6 

f(P,P,t) = /(arctan-,t), 
V 

and its regularity is now the one of the function (6,p) i— »• arctan ^. Therefore, 
the controlled dynamics of a test-particle is given by Hamiltonian : 

H c (e,p,t) = H lp (e,p,t) + f(9,p,t). 
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Figure 4. Poincare sections of a test-particle of Hamiltonian Hi p (9,p,t) (left) and 
H lp (8,p,t) + f(9,p,t) (right). 

In the (9,p) variables (see FigQJ, the control is successful in reconstructing 
some invariant tori around the macro-particle. In other words, it enlarges this 
regular structure. 



Conclusion 

In this paper, we considered a simplified mean-field approach to investigate 
the saturated dynamics of a Single Pass FEL. In particular, we showed that 
the size of the macro-particle can be increased by adding a small pertubation 
to the system, thus resulting in a low cost correction in term of energy. The 
main idea is to build invariant tori localized at specific positions : this method 
is utterly general and could be succesfully used to adjust the size of the macro- 
particle, thus possibly enhancing the bunching factor. Analogously, by limiting 
the portion of phase-space spanned by the macro-particle, one could aim at 
stabilizing the laser signal. Two future lines of investigation are foreseen. First, 
we intend to implement the computed control term in the framework of the 
original iV-body self-consistent picture and explore possible beneficial effects 
for the evolution of the radiation field. Further, it is planned to apply the 
above technique to the case of the reduced Hamiltonian of [2113] . We point 
out that an experimental test of the control method on a modified Travelling 
Wave Tube has been done [10J in absence of self-consistency. 
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